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in Fully Frustrated XY Models 
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The unbinding of kink pairs on domain walls in the fully frustrated XY model (on square or 
triangular lattices) is shown to induce the vanishing of phase coupling across the walls. This forces 
the phase transition, associated with unbinding of vortex pairs, to take place at a lower temperature 
than the other phase transition, associated with proliferation of the Ising-type domain walls. The 
results are applicable for a description of superconducting junction arrays and wire networks in a 
perpendicular magnetic field, as well as of planar antiferromagnets with a triangular lattice. 
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A fully frustrated (FF) XY model can be defined by 
the Hamiltonian 



H = - J^cos(</jj 
(y) 



(1) 



where J > is the coupling constant, the fluctuating 
variables ipi are defined on the sites i of some regular two- 
dimensional lattice, and the summation is performed over 
the pairs of nearest neighbors {ij) on this lattice. The 
non-fluctuating (quenched) variables Atj = —Aji defined 
on lattice bonds have to satisfy the constraint ^ Aij = tt 
(where the summation is performed over the perimeter 
of a plaquette) on all plaquettes of the lattice. 

For two decades such models (on various lattices) have 
been extensively studied in relation with experiments on 
Josephson junction arrays [Q, in which tpi can be associ- 
ated with the phase of the superconducting order param- 
eter on the i-th superconducting grain, and Aij is related 
to the vector potential of a perpendicular magnetic field, 
whose magnitude corresponds to a half-integer number 
of superconducting flux quanta per lattice plaquette. A 
planar antiferromagnet with a triangular lattice also can 
be described by the Hamiltonian (|l]) (with Aij = ±7r). 

The ground states of the FF XY models on square 
1^ and triangular |^ lattices are characterized by the 
U{1) X Z2 degeneracy, which suggests the possibility 
of two different phase transitions. One of them (the 
Berezinskii-Kosterlitz-Thouless transition ^-^) can be 
associated with unbinding of vortex pairs and the other 
with proliferation of the Ising type domain walls. 

Teitel and Jayaprakash Q have proposed that the tem- 
perature Ty of the vortex pairs dissociation cannot be 
higher than the temperature Tdw of the phase transition 
associated with domain walls proliferation. The argu- 
ments supporting this conjecture have been put forward 
in Refs. [|-10| and are related to the presence on corners 
of domain walls of fractional vortices, which are expected 
to screen the interaction of integer vortices at T > Tdvv- 

The application of the Hubbard-Stratanovich transfor- 
mation PI to the FF XY model on a square lattice al- 



lows one to reduce it [Q to the system of two coupled 
unfrustrated XY models, which in the limit of strong 
coupling becomes equivalent to the so-called XY-lsing 
model: 



(y) 



(2) 



where Si = ±1 is the auxiliary Ising type variable. In 
this model, the coupling of the phase variables ipi across 
any domain wall is completely absent. Although this 
property appears as a direct consequence of taking (with- 
out any justification) the strong coupling limit, and such 
a description fails to take into account the existence of 
fractional vortices, it has been suggested that the 
Xy-Ising model may turn out to be a reasonable ap- 
proximation for investigation of the FF XY model. 

In the present Letter, we demonstrate that, in the FF 
XY model on square or triangular lattices, the phase 
transition on a single domail wall, which takes place at 
Tk < Ty and consists in dissociation of pairs of logarith- 
mically interacting kinks |Q , induces for T > Tk the loss 
of phase coupling across the wall. This indeed makes the 
behavior of the FF XY model analogous to that of the 
XY-ls'mg model. We also show that the suppression of 
the phase coupling between different Ising domains leads 
to Ty < Tdwj Sit least when the phase transition associ- 
ated with domain wall proliferation is a continuous one. 

The FF XY model on square (or triangular) lattice 
being one of the simplest examples of a system with non- 
perturbative coupling between continuous and discrete 
degrees of freedom, the results are of interest not only 
in relation to experimental realizations mentioned above, 
but in a more general context of two-dimensional statisti- 
cal mechanics. In particular, we discuss in the conclusion 
their consequences for the interplay between the rough- 
ening and the reconstruction transitions We do 
not consider here the FF XY models on honeycomb px| ] 
and dice lattices, which are characterized by much 
more developed discrete degeneracies. 
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In the ground state of the FF XY model on a 
square lattice, the gauge invariant phase differences 
6ij — ipj — ipi — Aij on all bonds are equal [when reduced 
to the interval (— tt, tt)] to ±7r/4 in such a way that sum- 
mation of 9ij over the perimeter of each plaquette gives 
TTCT, where tr = ±1 is called chirality. The plaquettes 
with positive and negative chiralities regularly alternate 
with each other, forming the checkerboard pattern 
The discrete twofold degeneracy of the ground state cor- 
responds to the change of the signs of all chiralities. 

A domain wall can be defined as a topological ex- 
citation separating two ground states which cannot be 
transformed into each other by a continuous rotation. 
Schematically, it can be represented as a line on a lat- 
tice, each link of which separates two plaquettes with 
the same chirality (Fig. 1). A domain wall is character- 
ized by a finite energy per unit length; therefore at low 
temperatures all domain walls which appear as thermal 
fluctuations form closed loops. 

If one considers an infinite straight domain wall and 
fixes the state (the values of tfi) on one side of the wall, 
the state on the other side of the wall cannot be arbi- 
trary and depends on the position and on the orientation 
of the wall . If the wall is displaced by one lattice 
constant, the values of ipi on the other side of the wall 
are changed by tt. 

The presence of a kink (of the minimal height) on a do- 
main wall [Fig. 1(a)] produces a mismatch of tt between 
the states which have to be obtained when crossing the 
left and the right parts of the wall. This discrepancy has 
to be taken care of by a fractional vortex with the topo- 
logical charge ±1/2 located on the kink. The energy of 
such simple kink is therefore logarithmically divergent. 
The kinks of the double (or, generally, even) height [Fig. 
1(b)] do not introduce any mismatch, and their energy is 
finite. 
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FIG. 1. A domain wall with (a) a simple kink; (b) a double 
kink. Pluses and minuses show the signs of chiralities. 

Let us consider an infinite domain wall, introduced, 
for example, by an appropriate choice of boundary con- 
ditions. At low temperatures, it should contain a finite 
concentration of free double kinks, but all simple kinks 
have to form neutral pairs. Therefore, although the fluc- 
tuations of the domain wall diverge, the symmetry with 
respect to its shift by one lattice constant is broken. 

With an increase of temperature, the phase transition 
in the one-dimensional logarithmic gas of simple kinks 
will lead to dissociation of neutral pairs and to the ap- 
pearance of a finite concentration of free simple kinks 



[ p^ . As follows from the renormalization group analysis 
of Ref. [|l^, this takes place when the prefactor of the 
logarithmic interaction of simple kinks is equal to 2T, 

that is, at 

Tk = Jr(rK) . (3) 

Here T{T) is the helicity modulus, the macroscopic pa- 
rameter describing the effective stiffness of the system 
with respect to continuous twist of (p. At zero temper- 
ature r(0) = J/\/2, whereas the unbinding of integer 
vortices "in the bulk" takes place at [|l9| 

Tv = f r(Tv) , (4) 

that is, above Tk. The numerical evidence for unbinding 
of kink pairs at Tk < Ty has been obtained by Lee et al. 
[0, who, however, mistook this phase transition for the 
roughening transition of domain walls. 

The phase transition associated with unbinding of 
pairs of simple kinks leads to a restoration of the symme- 
try between the odd and the even positions of the domain 
wall and also to the loss of the effective phase stiffness 
across the wall. Any attempt to create a phase gradient 
perpendicularly to the wall will be relaxed due to motion 
of free simple kinks along the wall in different directions 
(in accordance with the sign of the topological charge) 
under the action of Magnus force. The situation is anal- 
ogous to what happens in the bulk above Ty, when the 
presence of free vortices prevents the creation of any sta- 
tionary phase gradient (supercurrent). For T < Tk all 
simple kinks are bound in neutral pairs and their relative 
displacement requires a finite energy, which means that 
the phase stiffness across the wall remains finite. 

Analogously, a phase gradient parallel to the wall will 
not penetrate on the other side of the wall. Instead there 
will appear a difference in concentration of simple kinks 
with positive and negative topological charges, which will 
compensate for the difference in phase gradients on both 
sides of the wall. Although creation of such difference in 
concentrations requires some energy, this energy is pro- 
portional to the length of the wall, whereas penetration 
of the phase gradient across the wall would require the 
additional energy which is proportional to the total area 
of the domain on the other side of the wall. The same 
happens on grain boundaries in crystals, where the dif- 
ference in orientations is taken care of by a sequence of 
dislocations of the same sign. 

Note that both mechanisms work only at length scales 
which are large in comparison with the inverse linear 
concentration of free simple kinks. Nonetheless, their 
existence implies that at large length scales the FF 
XY model indeed can be approximated by the XY- 
Ising model as proposed in Refs. p^ , p^ . Naturally, at 
T < Tdw such equivalence works only in a small vicinity 
of Tdw J in which the correlation length, defined in terms 
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of (J, is much larger than the typical distance between 
free simple kinks on a domain wall. 

The same conclusions are valid in the case of the an- 
tiferromagnetic XY model on a triangular lattice. The 
ground state of this model consists of three sublattices, 
the values of (fi in which differ by ±27r/3, and also is char- 
acterized by the C/(l) x Z2 degeneracy If the ground 
state on one side of a straight domain wall is fixed, the 
state on the other side of the wall can be obtained by a 
permutation of values of ifi on any two sublattices and 
subsequent rotation of all variables by tt . 

The three available options correspond to three posi- 
tions of the wall and differ from each other by global 
rotation by ±27r/3. Therefore, the simple kinks separat- 
ing the straight parts of a domain wall have to behave 
as fractional vortices with topological charges ±1/3. Ac- 
cordingly, the phase transition associated with kink pairs 
unbinding on an isolated infinite domain wall takes place 
at = (7r/9)r(T^), which is again below Ty. As in 
the case of a square lattice this phase transition leads 
to the loss of phase coupling across a domain wall and 
makes the behavior of the system analogous to that of 
the XF-Ising model. 

The helicity modulus F can be defined, in particu- 
lar, through the response of a system to application of 
the specially chosen boundary conditions (see, for exam- 
ple, Ref . ) . The important property of the XF-Ising 
model is that, as soon as there is at least one domain 
wall crossing the whole system, the variables ip at op- 
posite boundaries are completely decoupled from each 
other, which means F = 0. In the thermodynamic limit 
this takes place at any temperature higher than Tdw- 

At T < Tbw an externally imposed twist of ip in any 
typical configuration of Si can be carried only by the 
largest (infinite) cluster formed by the sites with the same 
sign of Si [note that the variable Si of the XY-ls'mg model 
(^ should be associated with the staggered chirality of 
the FF XY model and not with the chirality itself]. All 
other clusters have finite size and therefore are insensitive 
to boundary conditions. 

In two dimensions (in contrast to three), the point of 
the phase transition of the Ising model coincides with 
the percolation transition in the system of spin clus- 
ters so the density of the infinite cluster decays on 
approaching Tdw from below as p{T) cx £,p^'^- Here 
Ad = 2 — d = 5/96 |22| is the deviation of the fractal 
dimension d of the infinite cluster (at T = Tbw) from 
its Euclidean dimension and £,p{T) cx (Ibw ~ T)^'^ is 
the percolation length, the temperature dependence of 
which in the Ising model is described by the same expo- 
nent 1^=1 H^] as that of the correlation length ^. 

Therefore the bare (i.e., not reduced by the fluctua- 
tions of if) helicity modulus Fo(T) on approaching Tdw 
has to decrease algebraically: Fo(T) cx (Tdw — T)*, at 
least as fast as p{T) (actually much faster), which can 



be shown with the help of the variational calculation. 
The vortex pairs dissociation takes place as soon as the 
(renormalized) value of F is reduced to (2/7r)T, that is, 
below Tdw- For Tdw < t|°^ [where T^°^ w (7r/2)F(0) 
is the naive estimate for Ty], one therefore can expect 

Tdw^Tv cx [Tdw/Tv^V^* • "^^^ transitions can hap- 
pen simultaneously only if they occur as the first-order 
phase transition with larger than universal jump in F. 

According to our results the analogous behaviour can 
be expected in the FF XY models. However in that case 
the dependence of Tdw ~ on Tdw will be more com- 
plicated, since the effective reduction to the XY-lsmg 
model is developed only at the length scales which are 
large in comparison with typical distance between free 
simple kinks. The idea that Fo(T) is strongly suppressed 
on approaching Tdw is supported by a comparison of 
the results of numerical simulations of the same system 
at zero and full frustration [Q, which shows that in the 
latter case the drop of F with an increase of temperature 
for the same size of the system is more sharp. 

In conclusion, until now the use of the XF- Ising model 
for the description of the properties of the FF XY model 
on a square lattice |l^,|2^ could be considered a rather 
arbitrary procedure. Since application of the Hubbard- 
Stratanovich transformation [ ll] to the FF XY model on 
a triangular lattice is known to produce a coarse-grained 
Hamiltonian with a wrong symmetry of the ground 
state [U{1) X Z3 instead of U{1) x Z2], one always could 
have doubts if the application of the same transforma- 
tion on square lattice does not lead to the loss of some 
important properties of the original model. 

An argument in favor of such an approach has been put 
forward by Knops et al. 20 1, who have shown with the 



help of the numerical transfer matrix diagonalization that 
at T = Tdw the free energy of the 19- vertex version of the 
FF XY model with increasing the system size becomes 
insensitive to the boundary conditions inducing the twist 
of (fi. The authors of Ref. [|2^ have interpreted this as 
evidence for irrelevance at criticality of the coupling of 
ip across a domain wall. The same observation could be 
alternatively interpreted as evidence for Ty < Tdw- 

In the present Letter we have demonstrated that in the 
FF XY models on square and triangular lattices the loss 
of the phase coupling across domain walls is achieved al- 
ready at temperatures below Ty due to the presence on 
domain walls of free simple kinks. We also have shown 
that the dissociation of vortex pairs has to take place 
at Ty < Tdw, because, on approaching Tdw from be- 
low, the part of the system which reacts to external twist 
becomes more and more dilute [25 . This makes the sce- 
nario of a single phase transition with a novel critical 
behavior [9|j2^ impossible, at least for T > Tk. 

These conclusions are not dependent on the particular 
form of the interactions in the system (as soon as the 
degeneracy of the ground state remains the same), and 
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are applicable, for example, also when the interaction of 
further neighbors is taken into account pTf . They are in 
agreement with the results of the most recent numerical 
simulations of the FF XY models on square |2^,^ and 
triangular ||3^ lattices, as well as of the equivalent (half- 
integer) Coulomb gas and of the SOS-Ising model 
p2[ which is (partially) dual to the XY-lsing model. 

Recently Lee et al. have demonstrated that the 
generalized XF-Ising model is dual to the SOS-Ising 
model introduced by den Nijs for the coarse-grained 
description of the interplay between roughening and 
missing row reconstruction on a surface of a crystal with 
simple cubic lattice. The proposed phase diagram of this 
model can be found in Fig. 3 of Ref. For the case of 
A = (which corresponds to the XF-Ising model with a 
complete absence of phase coupling across domain walls) 
it contains two regions, one (at i? < 0) with separated 
and the other (at i? > 0) with coinciding phase transi- 
tions. Our analysis, as well as the results of the numerical 
simulations of Ref. |Q , implies that the two transitions 
should be separated for both signs of R. 

It seems worthwhile to mention that the conclusion on 
a larger than universal value of r(Tv) in the FF XY 
models p9|-|3l|] has been obtained with the help of the 
Weber-Minnhagen (WM) scaling analysis |3|], which is 
based on the same renormalization group equations |^] 
as the universal prediction (^ and, accordingly, does not 
even allow for a possibility of a nonuniversal value of 
r(rv). The results of Refs. [|9|-|l) therefore should be 
interpreted as evidence for deviation from the WM scal- 
ing. The reasons for such a deviation can be easily under- 
stood. It follows from our analysis that, even when the 
presence of vortex pairs (the only factor taken into ac- 
count in the framework of the WM analysis) is neglected, 
r in the vicinity of Tdw has to be strongly scale depen- 
dent, because the effects related with suppression of the 
effective stiffness across domain walls develop only with 
the increase of scale. Thus it is important not to confuse 
the two mechanisms for suppression of F. The method 
for plotting the data, which allows for checking in what 
interval of length scales the WM scaling really holds, has 
been recently suggested in Ref. |2^]. 
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